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Abstract 



Matrix interpretations generalize linear polynomial interpretations and 

have been proved useful in the implementation of tools for automatically 

proving termination of Term Rewriting Systems. In view of the successful 

use of rational coefficients in polynomial interpretations, we have recently 

generalized traditional matrix interpretations (using natural numbers in 

\^ \ the matrix entries) to incorporate real numbers. However, existing results 

Xyy • which formally prove that polynomials over the reals are more powerful 

rj* ' than polynomials over the naturals for proving termination of rewrite sys- 

O . terns failed to be extended to matrix interpretations. In this paper we 

get deeper into this problem. We show that, under some conditions, it is 
possible to transform a matrix interpretation over the rationals satisfying 
a set of symbolic constraints into a matrix interpretation over the naturals 
#-/-\ ' (using bigger matrices) which still satisfies the constraints. 

Keywords: Matrix and Polynomial Interpretations, Program Analysis, 
Termination. 

q ; 1 Introduction 

Constraint solving is an essential technique for the implementation of automatic 
verification systems. Many verification problems can be expressed as sets of sym- 
j^ ■ bolic constraints which have to be tested for satisfaction or even solved to give 

some explicit solution certifying their satisfaction. For instance, termination 
problems are often expressed as conjunctions of weak or strict symbolic con- 
straints like e >z e' or e >- e! between expressions e and e' coming from (parts 
of) the programs [6]. Automatic termination tools have to check these con- 
straints and eventually provide an appropriate certificate. A standard approach 
is using algebraic interpretations which translate the symbolic constraints into 
some kind of arithmetic constraints. 

Example 1 Consider the following Term Rewriting System (TRS) TZ J7J 0/: 

/(/(*)) -> /(s(/P0)) (i) 

/(<?(/(*))) -» X (2) 



♦Partially supported by EU (FEDER) and MICINN grant TIN 2007-68093-C02-02. The 
final publication of this paper is available at www.springerlink.com 



A proof of termination with dependency pairs can be easily obtained as follows 
JSJ'.' Consider the following rules (called dependency pairs,) associated to 1Z: 

F(f(X)) -► F(g(f(X))) (3) 

F(f(X)) -> F(X) (4) 

The following polynomial intepretation with rational coefficients 

[f](x) = 2x + 2 [g}(x) = ix+i [F](aO = z 

can oe used to prove termination of 1Z by showing that [I] > [r] for the rules 
(Q), (0) (where [I] and [r] are the interpretations of terms I and r which is ob- 
tained by structural induction), and [u] > [v] for the dependency pairs (0), CI. 

A recent and fruitful approach is using matrix interpretations [TO] , where the 
/c-ary symbols / are given parametric matrix functions [/], e.g., F\Xi + • • • + 
Fk%k + Fq, where the i^'s are (square) matrices of some fixed dimension n and 
Fq is an n-tuple. The variables Xi, . . . ,Xk are intended to range on n-tuples as 
well. In [10] . only natural numbers are used both in matrices and n-tuples. 



Example 2 The following matrix interpretation over the naturals 
[/](*) = (] [)*+( I) Iff](«) = 

[F](x) = 



1 1 
1 

can also be used for proving termination oflZ in Example^ Both interpretations 
have been automatically obtained by using MU-TERM J7] 



In [Hd], Endrullis et al.'s framework was extended to matrices containing real 
numbers in the entries. The adaptation was motivated by a number of recent 
theoretical works and experimental evaluations showing that polynomial inter- 
pretations over the rationals can be advantageously used instead of polynomial 
interpretations over the naturals [S] [TTJ [T3] . In [TTJ , syntactic conditions ensur- 
ing that, when dealing with linear polynomial interpretations, real coefficients 
must be used for addressing the corresponding termination problem, were given 
for the first time. The extension of these results to matrix interpretations over 
the reals failed [U[2]. 

Thus, the following question arises: Are rational numbers somehow unnec- 
essary when dealing with matrix interpretations? By examining the proofs of 
termination contributed to the International Competition on Termination^ by 
tools like Jambox [5], which makes extensive use of matrix interpretations over 
the naturals, and comparing them to the corresponding ones generated by tools 
like mu-term, which emphasizes the use of polynomials over the rationals [15] . 
one may notice that, in many cases, proofs of termination with polynomials over 
the rationals somehow correspond to proofs using matrix intepretations whose 
matrices have specific shapes. Examples [TJ and [5] illustrate this connection, 
which we substantiate in this paper: the bigger are the values of the (possibly 
rational) coefficients in the polynomial interpretation the more non-null entries 
are in the corresponding matrix coefficients. In this paper we investigate this 



1 See http : //termcomp . uibk . ac . at/termcomp/ 



phenomenon. In Section [5] we develop a notion of numeric matrix representa- 
tion which permits a representation of a natural number as a matrix of (smaller) 
natural numbers. This representation preserves the usual arithmetics of natu- 
ral numbers (addition and product). Therefore, we can think of such matrices 
as having a value; the value of the number from which they were obtained. 
Then, in Section [3] we show how to extend this process to transform a matrix 
of natural numbers into a (bigger) matrix of (smaller) natural numbers. As a 
consequence, we prove that every matrix of natural numbers can be represented 
as a bit matrix with entries in {0, 1} which still preserves its 'value' and arith- 
metic behavior. In Section [H we address the problem of representing arbitrary 
rational numbers as matrices of integer numbers. We argue that this is possible 
only for finite subsets of rational numbers. We investigate the use of nilpotent 
matrices (i.e., square matrices which become null after a finite number of self- 
products) as suitable devices for achieving this. In Section [5j we introduce a 
new generalization of matrix interpretations, which we call block-based matrix 
interpretations. Essentially, we view a matrix as structured into blocks of (sums 
of) constant or scalar matrices. In Section|51 we investigate how the satisfaction 
of different kind of constraints is preserved under the matrix transformations 
investigated in the previous sections. As a consequence of our results, we prove 
that TRSs which can be proved terminating by using matrix interpretations 
over the naturals can also be proved terminating by using a matrix interpreta- 
tion based on bit matrices. Furthermore, we show that, under some conditions, 
proofs of termination which are carried out by polynomial or matrix intepre- 
tations over the rationals can also be obtained by using matrix interpretations 
over the naturals (like in Example [2]) . Section [7] summarizes our contribution 
and concludes. 

2 Numbers as matrices 

In the following, we use the standard notations and terminology for matrices 
[H OH HE]- Given p,q£ N >0 and a set of numbers N (usually R, <Q>, or N), we 
write A £ N pxq to say that A is a matrix of p rows and q columns with entries 
Aij (E N (apx q- matrix for short). If p = q, then A is a square matrix. 

We investigate the representation of (rational) numbers as matrices of integer 
numbers satisfying some structural properties. Of course, we want to ensure that 
the representation preserves (part of) the algebraic structure of the considered 
numeric domain, in such a way that, for instance, the arithmetic operations and 
orderings among numbers (of the considered kind) can be implemented by using 
matrix operations and orderings. 

We view a rational number as a product p- for an integer number p € Z and 
a positive natural number q G N>o- First, we formalize a generic framework 
for representing real numbers as matrices: mappings p : R — >■ K mx " and p : 
jgmxn _^ jg> provide a representation of real numbers as matrices and vice versa. 

Definition 1 (Numeric matrix representation) Let N C R be a subset of 
real numbers, p,q £ N>o, and M C W xq be a subset of matrices. A numeric 
matrix representation is a pair (p, p) of mappings p : N — > M (called a numeric 
representation^ and p : M —> N (called a matrix valuationj such that po p = id. 

Let lpxq be the p x g-matrix all whose entries contain 1. A matrix C = cl pX q 
for some c € R (i.e., whose entries are settled to c) is called a constant matrix. 



The identity (square) matrix of size n is denoted /„. A matrix S = cl n for 
some c E K is called a scalar matrix. We consider the following numeric matrix 
representation. 

Definition 2 Let m,p,qE N >0 and AE«. pxq . We let 

1. p^ q {x) — m^-lpxq, i.e., each real number x is mapped to a constant 



pq 



p x q-matrix with entries 2i£, 

f -i pq 

v p y^ 9 A- ■ 
2. p m (A) = — — — — — —, i.e., each matrix A is mapped to the number which 

is obtained by adding all entries in A and then dividing this number by m. 

In the following, we prove some properties of p m which are used below. 

Proposition 1 Let m,p,q E N >0 and A, B e W xq . Then, p m (A + B) = 
p m (A) + p m (B) and p m (aA) = ap m (A) for all a E K. 

Proposition 2 Let p,q,r e N >0 , A 6 W xq and B e W xr . If B (resp. A) 
is scalar and q < r (resp. q < p), or B (resp. A) is a constant matrix, then 
Pp (AB)=p p (A) Pq (B). 

Propositions Q] and [5] entail the following. 

Corollary 1 Let p,q,r E N >0 , A e W xq and B E K« xr . // B (resp. A) 
is an additive combination of scalar and constant matrices, then p p (AB) — 
P P (A)p q (B). 

Corollary 2 Let n E N>o and A, B be n-square matrices. If B is an additive 
combination of scalar or constant matrices, then p n {AB) = p n (A)p n (B) = 

Pn(BA). 

If we consider only n-square matrices for representations, then p' n {x) = xl n 
could also be used with p n as a numeric matrix representation. 

Remark 1 (Use of vectors) Since vectors v can be seen as special matrices 
v E W ixl of n rows and a single column, the numeric matrix representation in 
Definition^ can also be used to represent real numbers as vectors. In particular, 
we get p™ xl (x) = xl n and p n (v) = ^" i = 1 "' . 

3 Transforming matrices of numbers 

We can extend any numeric representation p : R — > W xq to a mapping p : 
j^mxn _^ ^mpxnq f rom m x n-matrices A into mp x ng-matrices p{A) by just 
replacing the numeric entries Aij in A by the corresponding matrices p(Aij), 
i.e., p{A) = (^(Aj))l^i fci • The new matrix can be viewed as a block matrix 
whose blocks are p(Aij) for 1 < i < m and 1 < j < n. 

Example 3 We can transform the matrix I I by using pf^ q in Definition 

fJl' with m = p = q = 3, we obtain: 

3x3,/ 3 V _ / pl x3 (3) pl x3 (0) \_f l 3x3 3x3 



^ (i 3 J'-{pl x3 (0) pl x3 (3) J-\0 3x3 l 3 x. 



The interesting feature of the matrix obtained in Example [3] is that it is a bit 
matrix whose entries are eithei|j or 1. Note that this is due to the use of /ig X3 
which permits a representation of '3' as a constant matrix 13x3 with l's only. 
For the numeric matrix representation in Definition [51 we have the following. 

Proposition 3 Letm,n,p,q,r,s G N> and A G W xs . Then, p mn (/j,f£ q (A)) = 

Pn(A). 

Theorem 1 Let m,n,p,q,r,s,t,u G N >0 - If A,B G R rxs , then, p m (A) + 
p m (B) = Pmn (pP l x "(A)+^ q (B)) and Pm {aA) = ap mn {pP n Xq {A)) for alia € R. 
IfAe R sxt and B G R txu , then p m {AB) = p m n(p% xq (A) fj^ r (B)) . 

Propositions [1] and [5] entail the following. 

Corollary 3 Let m,n,p,q,r,s G N> , A G R mxm and B G W nxs . If A 
is an additive combination of scalar and constant matrices, then p m (AB) = 
p m M xq (A)p? xr (B)) - Pmn « /q {A)) Pmn {pT r {B)). 

3.1 Representing integer numbers as matrices 

In the following, we use /i^ x<? in Definition [5] as a basis for the following repre- 
sentation mapping for integer numbers. 

Definition 3 (Representing integer numbers as matrices) Let n G N be 

such that n > 1 and fi n be given as follows: for all x G Z ; 

— l„xn if n divides x 



1 \ f flnxn if n divide 

Pn(x) = < n T ±1 

{ xl n otherwise 



We also define v n (x) — xl n to represent a number x as a n-dimensional vector. 

( 2 2 \ 
Example 4 The matrix p,2 (4) = I I represents 4 according to Definition 

m 

Note that, for all n G N, (/ii ra , /?„) (with p n as in Definition [5]) is a numeric matrix 
representation for integer numbers. We obtain the following: 

Proposition 4 (Decrease of natural entries in matrices) Every matrix A G 
fqpxq suc j l that n — max(A) > 1 can be represented as a matrix A' G f$ n P xn i 
such that, for all m G N>o, p m (A) = p m n{A') and max(A) > max(A'). 

Obviously, Proposition 2] entails the following. 

Corollary 4 (Natural matrices as bit matrices) Every matrix A G N pxq 
can be represented as a bit matrix A' G {0, 1}™p x ™9 for some n G N>o and for 

all m G N>o, p m (A) = p mn (A'). 



2 Matrices with entries in {0, 1} are called (0, l)-matrices in 1181 Section 8.2]. 



4 Representation of rational numbers below 1 

In this section, we investigate matrix representations (over the naturals) which 
can be used to deal with rational numbers - for some q G N>o- A ^-square 
matrix Ax which is almost null except for a single entry of value 1 can be used 

to represent - because p q (Ai) = -. By Corollary [21 for A p = pl q we get 
p q (ApAi) = p q (AiA p ) — p q (A p )p q (Ai) — 2, Therefore, we can represent a 
rational number ^ asa g-square matrix with a single entry of value p. However, 

we have to change the size of the matrix if a different number "P-, with q ^ q' is 
considered. 

Remark 2 Note that there is no generic representation of all rational num- 
bers by using matrices over the naturals of a given dimension which is able to 
represent their values and appropriate comparisons among them. Such a rep- 
resentation should be able to represent a decreasing sequence 1 > \ > | > 
• • • > — > • • • by means of matrices Ax G N pX9 for all n G N>o satis- 
fying p m (Aj_) > p m {A_L_) for all n G N>o- Equivalently, we should have 
mp m (A±) > mp m (A^^) for all n G N>o- Since mp rn {A±),mp m (A_i_) G N. 
this would imply the existence of an infinite decreasing sequence of natural num- 
bers, which is not possible. 

Furthermore, the product of rational numbers - and — represented as the q- 
square matrices indicated above is not preserved by the matrix product. Thus, 
in the following, we look for better representations. 

4.1 Use of nilpotent matrices 

Nilpotent matrices are n-square matrices B satisfying B k — 0„ xn for some 
positive integer k G N>o [IB Section 4.1] (which is called the degree of nilpotency 
[HI page 12] or the index of nilpotency of B [TTJ page 396]). The degree of 
nilpotency k of a n-square matrix A is bounded by n: k < n d7, Exercise 7.7.1]. 
Given n, the following n-square matrix (called a Jordan block JT7J Page 579]): 

/0 1 \ 



'■• 1 

V o/ 

(i.e., there are n — 1 ones in the superdiagonal of J n and all other entries in J n 
are zeroes) is nilpotent of degree n, i.e., J™ = nX n- 

Write J n — [0, Z\, . . . , Z n _i], where Zi is an almost null vector containing 
a single 1 in the i-th entry. For instance, Z\ = (1,0, .. . ,0) T . Then, it is not 
difficult to see that J% is obtained by introducing p — 1 columns of zeros from 
the left side of the matrix and shifting columns Z; to the right just throwing 
those which exceed the n-th position, i.c, J? = [0, . . . , 0, Z±, . . . , Z„_ p ]. 

Remark 3 J n Is also known as a shift matrix because, for an arbitrary matrix 
A, J n A is obtained by shifting the rows of A upwards by one position and in- 
troducing a new bottom row of zeroes. Similarly AJ n shifts the columns of A to 
the right and introduces a new leftmost column of zeroes. 



Note that, for all p € N, 

I II T) "^ Tl 

The following result is obvious. 

Proposition 5 Lei m,n E N>o a«rf < p < q < n. Then, p m (Jn) = 
Pm((JZT) > p m {{Jl) q ) = Pm {J q n ). For all A E R nxr , p m (J%A) > p m {J q A). 
For all A E W* n , p m {AJP) > Pm (AJ%). 

In general it is not true that p m (J p A) = p m (AJ p ) for square matrices A. Due 
to Proposition [5] and Corollary [TJ for additive combinations A of constant and 
scalar matrices we have the following: 

Corollary 5 Let m,n E N> awd < p < q < n. If A e W lXr (resp. 
A G W xn ) is an additive combination of constant and scalar matrices, then 
Pm {.P n A) = p m (JP)p m (A) > p m {J q n ) Pm {A) = p m (J q A) (resp. p m (AJ%) = 
Pm (A) Pm (JP) > p m (A) Pm (J*) = Pm {AJD). 

Thus, the following property of rational numbers r = — for some n € N>o is 
simulated by the representation of integer numbers as (additive combinations 
of) constant or scalar matrices, and rational numbers — as powers of Jordan 
blocks: for all n £ N and positive rational number r, < r < 1, we have 
n > nr > nr 2 > ■■ ■ . 

Example 5 We can use Ji to represent ^: p2^J , i) = ^- However, ,]%, which is 
expected to correspond to j does not fit this expectation: P2(Ji) — /^(Cb) = 0. 

Theorem 2 If A\, ■ ■ ■ , An are n-square matrices such that, for all i, 1 < i < 

N , either 

1. Ai is scalar or constant, or 

2. Ai — J Pi for some p b E N and then both Ai-\ (if i > 1) and Ai + \ (if 
i < N ) are constant matrices, 

then PnflliLl A i) = n^Il Pn(Ai). 

As remarked above, it is not possible to use matrices of natural numbers of a 
fixed size n to represent all fractions - for q E N>o- Instead, we will consider 
the problem of representing finite subsets Q C {i | q E N>o} by using n-square 
(nilpotent) matrices in such a way that the following property is fulfilled by the 
representation (p,, p): for all x,y E Q such that xy E Q, p(p(x)p(y)) — xy — 
p{p{x))p(p{y)), i.e., the number p(p,{x)p(y)) which corresponds to the matrix 
product p(x)p.(y) of matrices p(x) and p{y) representing x E Q and y E Q is 
exactly xy E Q. In Section [6] we discuss how to take benefit from this. 

The dimension n of the matrices involved in the representation of Q is de- 
termined by the least element in Q. For instance, we can fix n such that — is 
the least number in Q. Then, an obvious representative for — is J^~ l because 
/? n (J™ _1 ) = — (see Example[S]). However, the feasibility of this simple approach 
depends on the other values in Q. 



Example 6 Let <2 = {4, 4}. If we fix J| to be the representation of \, then the 
representation of \ should be J\ (because P4(J|) = \). However, p4((J|) 2 ) is 
not j as one could expect; instead, /^((Jf) 2 ) = 0. The following block matrices 
of size 4 whose blocks are combinations of (transposed) Jordan blocks can be 
used to represent Q as required: 

Qi = ( J2 ^ ") Qi = ( ° 2x2 J2J z 

2 \ 02x2 02x2 / 4 \ 02x2 02x2 

Note that p4Qi) = i±i = ±, p 4 (Q i) = ^ = \, and Q1Q1 = Qi 
Example 7 Let Q = {|, j, i}. The following matrices of size 8: 

J& J4 \ p. f J4 04x4 

04x4 04x4 / T V 04x4 04 X 4 



Qi = 



J4 04x4 

04x4 04x4 



can be used to represent Q. Note that ps(Qi) = g - = ?' P&(Q±) ~ ~ s = !-• 



2 i P»\>«i/- 8 ~ 4> 

__ Di Di — Di r,T,rin,n, —DiDi — 

8 8' 



PsiQi) = 8^ = g, Q1Q1 = Qi and Q1Q1 = Q1Q1 = Qi, as required 



Example 8 Lei Q = {5, i, A}. TTie following (block) matrices of size 6 can &e 
used to represent Q: 

h (JIV \ Ql = ( J i MJi) T 

03x3 O3X3 / 3 \ O3X3 03 x 3 



Ql 



03x3 -^(-^V 
03x3 3X 3 



Note that p 6 (Qi) = 2±1 = I ; ^(Qi) = 1±1 = I, an d ^(Qi) = 2±i = I. 
Again, it is not difficult to see that Q1Q1 — Q1Q1 — Q±. Note, however, that 
if we add j to Q, then we could not use these matrices for representing j. In 

particular, Q\ = Qi, i.e., pe{Q 2 i) 7^ j as should be the case. 

2 3 2 

5 Matrix interpretations revisited 

As remarked above, termination problems in term rewriting are usually trans- 
lated into conjunctions of weak or strict symbolic constraints like s >z t or s >- 1 
between terms s and t coming from (parts of) the TRS. In order to check the 
satisfaction of these constraints, we need to use term (quasi-) ordcrings. Such 
term ordcrings can be obtained by giving appropriate interpretations to the 
function symbols of a signature. Given a signature T ', an F-algcbra is a pair 
A = (A, .Fa), where A is a set and Fa is a set of mappings fj\_ : A k — > A for 
each /£ J where k = ar(f). For a given valuation mapping a : X — > A, the 
evaluation mapping [a] : T(J-,X) — > A is inductively defined by [a](x) — a{x) 
if x e X and [a]{f(t 1 ,...,t k )) = f A ([a]{ti), . . . ,[a]{t k )) for x <= X, f e F, 
t\, . . . , tk G 7~(F, X). Given a term i with Var(t) — {cci, . . . , x„}, we write [£] 
to denote the function F t : A™ — > A given by F t (ai, . . . ,a n ) = [a( ai a n )](£) 
for each tuple (01, . . . , a n ) € A™, where av ai .... i a„)( 2; i) = a i f° r 1 < « < n. We 
can define a stable quasi-ordering >; on terms given by t >^ s if and only if 



[a](t) ^a [a](s), for all a : X — > A, where ^a is a quasi-ordering on A. We can 
define a stable strict ordering Zl on terms by t Zl s if [a](t) >~a H(s), for all 
a : X — > A, where ^a is a strict ordering on A. 

A matrix interpretation for a fc-ary symbol / is a linear expression F\X\ + 
• • • + -FfciEfc + Fq where the F%, . . . , Ff. are (square) matrices of n x n natural 
numbers and the variables x\,...,Xk (and also the constant term Fq) are n- 
tuples of natural numbers [TU] ■ An expression like Fx, where F is an n-square 
matrix and x is an n-tuple of numbers, is interpreted as the usual matrix- vector 
product, i.e., the i-th component yi of y = Fx is yi = Xo=i ^ijXj. Matrices 
and vectors are compared by using an entrywise ordering: for A,B€ W xq 1 we 
write A > B if Ay > By for all 1 < i < p, 1 < j < q [H Chapter 15]. We also 
write A > B if A > B and An > B\\ [TU]. Note that this also describes how to 
compare tuples of numbers. In [TJ[2], Endrullis et al.'s approach was extended 
to matrix interpretations with real numbers in matrix entries. 

Here, we generalize the notion of matrix interpretation in the following ways: 

1. We consider a domain T n ^(N) of block-based tuples over N , i.e., n-tuples 
consisting of /3 = ? tuples of size b £ N>o which are constant tuples clt 
for some number c € N. If we take 6=1 and N = N or N = Rq, then we 
are in the original approaches [TUJ and [TJH], respectively. Note, however, 
that given n £ N>o, only divisors 6 of n can be used to establish blocks 
within n-tuples of numbers. 

2. Matrices F in matrix expressions interpreting symbols / £ T will be block 

/ Fn ■■■ F 1/3 \ 

matrices : • . : such that i 7 ^ = Cy + Sy is a sum of a 6- 

V Fpi ■■■ F^ J 
square constant matrix Cy = Cijlbxb, where cy £ iV, and a 6-square 
scalar matrix SV,- = Sijlf,, where sy £ AT, for all 1 < i,j < f3. This 
is necessary for soundness of the obtained algebraic interpretation: the 
product of one of such matrices by a block-based tuple as above produces 
a block-based tuple as above, i.e., Fv £ T n ^{N) for all v £ T ny b{N). 
Furthermore, such matrices are closed under addition and matrix product. 
This is essential to obtain matrices of this kind during the interpretation 
of the terms, where nested symbols yield products and sums of matrices 
after interpreting them. Again, if b = 1, we are in the usual case for matrix 
interpretations. 

3. Given a matrix valuation p, and matrices A, B, we let A > p B if p{A) >n 
p(B), where >jv is an ordering over N. Given <5 > 0, the following ordering 
over (real) numbers is used 115) : x >s y if x — y > S. If 8 = 1 and x,y 
are natural numbers, we obtain the usual well-founded ordering among 
natural numbers >n- Now, the following (strict and well-founded) ordering 
>p,8 (or just > p or even > if it is clear from the context) on n-tuples of 
nonnegative numbers is considered: x > p ,s if if p{x) >s p{y)- Clearly, 
(T! n ^(N), > p ,i) and (T n ^(S.o), >p,s) are well-founded orderings. 

The previous orderings do not take into account the block structure of 
matrices or tuples. The following one does: for matrices A, B with blocks 
Aij and -By for 1 < i,j < (3, we write A > b B if Aij > p Bij for all 
1 < i,j < P- Similarly, A > b pS B if An > Pt $ B lx and A tj > p By for all 
1 < z , j < /3 . These definitions are adapted to tuples in the obvious way. 



Definition 4 (Block-based matrix interpretation) Let T be a signature, 
n, b G N>o be such that b divides n, and j3 — 7 j-. An (n,b)-block-based matrix 
interpretation is an J- -algebra A = (A, J- a) such that 

1. A = T n . b {N) for N = N or N = R , and 

2. J- a consists of matrix functions [f](xi, . . . , Xu) — F\X\ + • • • + -Ffc-Xfc + Fq 
which are closed in A and where, for all f G T , 1 < % < k, 

(a) Fi is an n-square matrix of /3 x /3 blocks of b- square matrices Cji+Sjz 
such that Cje is a constant matrix, and Sjg is a scalar matrix for all 
1 < j < 13 andl<£< [3. 

(b) Fq — (cf • • ■ Ca) T consists of (3 b-tuples Cj — Cjlb f or some Cj G N. 

We make the orderings which are going to be used explicit by adding them to the 

P J, (A,F A ,> b J, 



pair {A, Fa), thus specifying an ordered algebra: {A,Fat>p„), (A,J 7 A,> b 



etc. 



6 Solving and transforming matrix constraints 

In the following, we consider two kinds of constraint solving problems which are 
relevant in proofs of termination. 

6.1 Testing universally quantified symbolic constraints 

Proofs of termination in term rewriting involve solving weak or strict symbolic 
constraints s >z t or s Zl t between terms s and t coming from (parts of) the 
rules of the TRS where the variables in s and t are universally quantified in 
the corresponding constraint. Here, y and Zl are (quasi-)orderings on terms 
satisfying appropriate conditions [3) [6j |8] . 

Example 9 The following symbolic constraints must be checked in order to 
guarantee termination of the TRS 1Z in Example \^ 

VX(/(/(X)) h f(g(f(X)))) (5) 

VX(f(g(f(X))) t X) (6) 

VX(F(f(X)) zi F(g(f(X)))) (7) 

VX(F(f(X)) Z F(Xj) (8) 

Here, variables X range on terms in T(T,X\ >^ and Z are intended to be in- 
terpreted as a monotonic and stable quasiordering on terms, and a well-founded 
and stable ordering on terms, respectively. 

In the so-called interpretation method for checking symbolic constraints as the 
ones in Example^ we use appropriate ordered J- -algebras to generate the neces- 
sary orderings (see Section [S]). In our setting, we are interested in investigating 
the use of matrix algebras. The following result shows how to decrease the 
value of some of the entries in a block-based matrix algebra over the naturals 
to obtain an equivalent block-based matrix algebra which uses bigger matrices. 



Theorem 3 Let T be a signature and c be a symbolic constraint \fx(s cxi t) 

>b ) 
-Pb> 



for terms s,t G T(T,X). Let A = (T!„.b(N), J\4, >^ ) be an (n,b) -block-based 
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matrix interpretation over the naturals and N be the maximum of all entries 
occurring in any matrix or vector in A. Let B = (B, J-g, > b J? ) be an (Nn, Nb)- 
block-based matrix interpretation where B — 2jv>i,jv&(N) and for all f G T , 

[/]s(*i, •■•,Xk) = m( F i)xi H \-m(Fk)xk + vn(F ) iff [/]^(*i,.. . ,Xk) = 

FiX\ + • • • + FkXk + Fq. Then, A satisfies c if and only if B satisfies c. 

Example 10 Consider the following TRSs fl(A Example 3}: 

K:f(a,g(y),z) -> f(a,y,g(y)) (9) 

f(b,g(y),z) -► f(a,y,z) (10) 

a -> b (11) 

S:f(x,y,z) -»• f(x,y,g(z)) (12) 

In order to prove termination of 1Z relative to S (written SN(R/S)), we have 
to check that \/x(l >; r) holds for all rules I — >• r G 1Z U S. Endrullis et al. use 
the following matrix interpretation for that: 



[b] 



[f](x,y,z) = (^n) a; +(nn)y+(nn) z + 



[<?](*) = 



1 


\ 





) 


' 


\ 





) 


1 











1 





1 


1 



/ y ' V ' 



By using Theorem^ we conclude that the following block-based matrix inter- 
pretation with bit matrices of dimension 4 can also do the work. 

fa 

0a 

,i / h 02x2 ^ . ( h 1-2x2 \ . / h 2X 2 "\ . ( O2 

W(*.»'*) - (02X2 02X2 X+ { 02X2 02X2 j^l 02X2 02x2J Z+ (o2 



[6] 



02 
02 



^2 2X 2 



ww = ^ /rj- + V i2 

6.2 Solving existentially quantified arithmetic constraints 

In many applications, when algebraic interpretations have to be generated rather 
than given by the user, it is usual to work with parametric interpretations. 

Example 11 For the symbols occurring in 1Z in Example^ we can consider 
the following linear parametric interpretation: 

[f](x) = fix + f [g](x) = gix + g \F](x) = F x x + F Q 

where fi,gi, and Fi are expected to be n-square matrices for some n G N>o (the 
case n = 1 corresponds to a linear polynomial intepretation) and fo,go,Fo are 
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n-tuples. The satisfaction of the following arithmetic constraint^ 



hh 


> 


figih 


/1/0 + fa 


> 


/1.91/0 


/iSi/i 


> 


1 


/lffi/o + fi9o + fa 


> 





Fifi 


> 


Figih 


F\fa + F 


> 


Figifo 


Fifi 


> 


Ft 


F\fa + F 


> 






(13) 
fi9o + fa (14) 

(15) 
(16) 
(17) 
- F m + F (18) 

(19) 
(20) 

is necessary to ensure that 1Z is terminating. 

By a parametric matrix interpretation we mean a matrix intepretation where 
the entries in matrices are not numbers but rather parametric coefficients, i.e., 
variables for which we have to provide a numeric value, depending on some 
(existential) constraints. 

In general, given a set of variables X, we consider here symbolic arithmetic 
constraints of the form s cxi t for cxi G {>,>}, where s is of the form X^=i s * 
for m s > and Si = su • ■ • Si ms t with m S] j > for all 1 < i < a s and Sy G X 
(t would have an analogous structure). Furthermore, we assume existential 
quantification over all variables occurring in s and t. Note in Example 1 1 1 1 that 
we use constants like and 1. They are handled as 'special' variables which will 
receive the intended interpretation. 

Consider a valuation r\ : X — > N for the variables in X as numbers in 
N . Here, wc consider N = N U Q for some finite subset of rational numbers 
Q Q Q— N satisfying some conditions. We are interested in representing numbers 
•q(x) as matrices /j,(r](x)) as discussed above. Note that we cannot represent 
arbitrary rational numbers by using matrices over the naturals (Remark [2]) . 
Still, when dealing with finite sets C of arithmetic restrictions, we can restrict the 
attention to those rational numbers which are required to check its satisfaction 
for a given valuation r\. This includes not only rational numbers 77 (x) which 
are assigned to x G X, but also those which could occur during the evaluation 
of an arithmetic expression due to products r](x)r](y) of rational numbers r)(x) 
and r/(y) which have been associated to variables x and y. The idea is that Q 
should contain such rational numbers. 

Definition 5 (Compatible domain of rational numbers) Let r\ : X — >• NU 

Q be a valuation for some QCQ-M and C be a set of arithmetic constraints. 
Given a multiplicative component s — Si ■ • • s m of an arithmetic expression in 
a constraint in C, let I s = {i±, . . . ,ifc} be the set of indices of variables in s 
whose valuation is a rational (and noninteger) number, i.e., for all i G I s , 
il(si) G Q — N. We say that Q is compatible with C and r\ if Y[i£jV( s i) £ Q 
for all multiplicative components s in C and J G I s such that J 5= 0. 

Example 12 The numeric valuation which corresponds to the poynomial in- 
terpretation in Example [7] is: 



77(1) 1 77(0) = 77(A) = 2 r,(/o) = 2 

\ v(go) = \ 



v(gi) = \ v(go) - h v(Fi) = i v(f ) = o 



3 By lack of space, we cannot explain how these constraints are obtained. Full details about 
this standard procedures can be found elsewhere, see, e.g., 7, 10, 15 16 . 
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The set Q = {i} is compatible with this valuation and with the constraints C in 
Examvle llll Consider C = CU{f±gifi > /i<?i/i3i}- Now, Q is not compatible 
with C and r\ because we have that r](gi)r](gi) = h ^ Q. If we add \ to Q, then 
we get compatibility with C . 

A valuation r\ is extended to terms and constraints by n(s X t) = rj(s) ixq T)(s), 
■q(si H h s m ) = r)(si) H h rj(s m ), and rj(xi ■ ■ ■ x m ) = j](x\) ■ ■ ■ rj(x m ). 

Remark 4 In general, we assume that + is commutative, but (as happens with 
the matrix product) we do not assume commutativity of the product in arithmetic 
expressions or their valuations. 

Now we have to extend p n in Definition [3] to deal with rational numbers in Q. 

Remark 5 In contrast to natural numbers, we have no systematic way to asso- 
ciate matrices to rational numbers yet. In Section R] we have investigated some 
partial solutions to this problem. In particular, Examples El and\^ show 
that the dimension n of the considered matrices and tuples heavily depend on 
the numbers in Q. On the other hand, these examples also provide useful en- 
codings for rational numbers which are frequently used in automatic proofs of 
termination with polynomial or matrix interpretations £31 1 1 1[ \16\j . 

Assume that /j, n has been extended to each x G Q in such a way that: for all 
x,y £ Q such that xy e Q, p n (Hn(x)n n (y)) =xy = pn(jtn(x))p n (pn(y))- 

In our setting, variables in X are interpreted not only as matrices but some of 
them as vectors. Assume that X$ C X must be interpreted in this way and that 
the constraints in C are consistent with this, i.e., whenever a variable xq G Xq 
occurs in a constraint c e C of the form s ix t, each multiplicative term in s and 
t must contain a single variable yo G Xq which must be at the end of the term. 

Example 13 For the constraints in Examvle Alll we have Xq = {/o, -Foj go, 0}. 

The vectorial (or n-tuple) representation of x £ N U Q by p, n is /j, n (x)l n . 

Example 14 The matrices over the naturals which correspond to r\ and Q in 
Examvle \lH\ is (with the encoding of ^ in Example[5j) are: 

til(fl)) = 12X2=(} }) tiv(fo)) = 2 1 2 = 

tivigi)) = J2 = f J tiv(9o)) = J2U = 

tiv(Fi)) = h=[\ ? ] tiv(Fo)) = 1 2 = 



1 

1 
1 



tivO-)) = h = ( I ? ) tiv(0)) = 1 2 = 




1 


1 


1 


1 


1 








1 



Accordingly, the matrix interpretation over the naturals which corresponds to 
the polynomial interpretation over the rationals in Example [7] is: 

[F}(x) 
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This matrix interpretation induced from the polynomial interpretation over the 
rationals in Example]]] can also be used to solve the constraints in the example. 

Our technique could be used to translate matrix interpretations over the ratio- 
nals into matrix interpretations over the naturals by just applying the previous 
translation to the entries of the matrix interpretation over the rationals instead 
to the coefficients of the polynomial interpretation. 

7 Conclusions 

We have investigated matrix representations of natural and rational numbers 
which can be used to simulate the arithmetics of natural and rational numbers, 
respectively. We have introduced the notion of numeric matrix representation 
(Definition [1]) which associates a number to a matrix and viceversa and proved 
that, by using some specific representations (Definitions [5] and [3]) , the arith- 
metic of natural numbers is preserved. Furthermore, we have proved that every 
matrix interpretation over the naturals has an associated bit matrix interpre- 
tation of (usually) bigger size of the same associated value (Corollary H]). We 
have investigated the representation of rational numbers by using matrices of 
natural numbers. We have proved that this problem has no general solution but 
we have found some suitable trade-offs for finite subsets of rational numbers by 
using nilpotent matrices consisting of Jordan blocks. Then we have introduced 
the notion of block-based matrix interpretation (Definition 0]) which generalizes 
existing approaches to matrix interpretations. We use it to transform matrix 
interpretations over the naturals into matrix interpretations over {0, 1}, and 
also to transform matrix interpretations over the rationals into matrix interpre- 
tations over the naturals. 

The question posed in the introduction: are rational numbers somehow un- 
necessary when dealing with matrix interpretations? could not be answered in 
full generality due to the lack of a general procedure for building matrix rep- 
resentations for arbitrary finite sets of rational numbers. Of course, this is a 
main topic of further research and we think that we have settled a good starting 
point in considering the use of combinations of Jordan blocks as in Examples 
El El El and M Nevertheless, our results suggest that the use of matrices over 
the naturals of big size can somehow play the role of rational numbers in inter- 
pretations of smaller size, and in particular in linear polynomial intepretations 
over the rationals. This does not mean that implementing polynomial or matrix 
intepretations over the rationals is not useful anymore and that natural num- 
bers should be used everywhere. In fact, working with matrix interpretations of 
big size is computationally expensive. Another interesting consequence of our 
analysis is the connection between dimension of matrices over the naturals and 
value of their entries via Proposition@]and Corollary@] Roughly speaking, these 
results can be interpreted by saying that bigger dimensions of matrices permit 
the use of smaller entries. In practice, most tools put strong numeric bounds 
to the coefficients or entries of the interpretations. Our results suggest that 
increasing such bounds could have a similar effect to increasing the size of the 
matrices. A more precise analysis about the trade-offs in design and efficiency 
which these considerations could lead to is also subject for future work. 
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